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A three dimensional Particle-In-Cell simulation describing the interaction of an intense laser beam 
with a plasma slab is presented. It is observed that the laser generated electron current decays into 
magnetically isolated filaments. The filaments grow in scale and magnitude by a merging process in 
the course of which the field topology changes. The opposite process also takes place occasionally. 
The laser driven charge and energy flows and the reconnecting magnetic field mutually interact. At 
the end of the merging process flows and fields are confined close to the laser irradiated surface of the 
plasma slab. Both decay rapidly in the bulk plasma. Due to the magnetic pressure in the filaments 
hollow density channels in the electron and ion densities are formed. The simulation reveals that 
charge flows in these channels can exceed the Alfven current. 
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■ Key issues of many applications concerning the interaction between lasers and matter at high intensities are laser 
Oh| absorption and charge and energy transport through a plasma. A typical application of intense laser-matter interaction 

where this is of relevance is Fast Ignition (FI) in Inertial Confinement Fusion (ICF) In the present paper we 

study issues of transport in a thick plasma slab with the help of Particle-In-Cell (PIC) simulations over a few hundred 

■ femtoseconds simulation time. Binary collisions are included. Details of the collisional model will be presented 
elsewhere. The density of the plasma slab is about 20 times over-critical. 

A reasonable definition of thickness has some ambiguity. We call a slab thick when its depth exceeds many times the 



^ [ skin lengths l s = c/lu p where w p is the plasma frequency. In addition, the thickness of the plasma slab must be larger 
• . than the penetration depth of the charge flow. The plasma slab investigated in this paper has a thickness which is 
large enough to exceed the latter sufficiently. Two dimensional (2D) transport simulations in the plane perpendicular 
r"| i to the laser direction have recently been reported B. In these simulations the slab thickness vanishes, the laser is 
Oh neglected and fast particles have been injected by hand. In addition, 2D geometry severely limits the available degrees 

of freedom for current transport and magnetic field evolution. 
y—{ [ In the present paper we investigate properties of laser-generated charge flows in a thick plasma slab in 3D where 
■ the fast particles are generated by laser irradiation. We show that magnetic filament merging [|| is an important 
process to form large filaments and magnetic field strengths. In 3D the process represents a simple form of magnetic 
reconnection in the sense that the field topology changes in the course of the latter. Laser generated charge flows 
found in the simulation exceed the Alfven limit . Furthermore, it is observed that the magnetic fields dominate 
i the properties of current transport. Both, electric and magnetic fields are localized close to the front surface of the 
t— I ! plasma slab as is the laser generated charge flow. 

The slab has a thickness of 6.0 /xm and a width of 4.0 /xm x 4.0 /xm. The box size is 4.0 /im x 4.0 /im x 8.0 /im. All 
fields depend on x, y, and z. The numerical grid has 152 x 152 x 800 cells. Lateral directions are periodic. Electrons 
. ' an< ^ ^ ons are P resen ted by 1.44- 10 s quasi-particles. The initial electron and ion temperatures are 10.0 keV and 1.0 keV 
respectively. The laser beam propagates in z and is linearly polarized along x. After a rise time of three optical cycles 

^ : 



the laser intensity is kept constant. The incident laser beam has a Gaussian envelope laterally with a width of 2.0/xm at 
full-width-half-maximum. The irradiancc is IX 2 — 5.0 ■ 10 19 Wcm~ 2 /xm 2 . The slab has a marginal initial deformation 
to enhance absorption The deformation is parameterized by z(x, y) = (5-cxp (—(a; — xq) 2 /r 2 + (y — yo) 2 /f 2 ) where 
8 = 0.4/im, xo = yo = 2.0 /xm and r = 2.5 /xm. The center of the slab is located at z = 4.2 /xm. The background 
plasma consist of protons. The initial electron plasma density in the simulation is ?x c = 3.33 • 10 22 cm~ 3 . The laser 
radiation is turned off after 200 fs. The simulation itself is stopped after 400 fs. The total simulation is carried out on 
a parallel computer using 361 compute nodes and consumes 29000 CPU hours in total. The coordinate system used 
in the simulations is right handed. Hence, the positive z-axis points out of the x/y-plane shown in the figures. In what 
follows wc will always use the positive z-axis for reference. 

Plots (a,b,c) of Fig. [I] show the plane z = 2.0 /xm of the cycle averaged magnetic field |B| defined in the figure 
caption. Time proceeds from plots (a,d) to (c,f). Filaments of different scale are observed. The white arrows indicate 
the direction of the magnetic field. It is seen that the magnetic field lines in the filaments are closed. As is seen from 
plots (a,b) the filaments attract each other. In case the attractive forces between them are strong enough magnetic 
field lines start to reconnect to form larger filaments with again closed field lines. The topology of the magnetic field 
changes. The inverse process takes place occasionally. The largest filament scale in our simulation is shown in plot 
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(c) of the figure. The peak magnetic field strength in this filament is about 45000 T. Plots (d,e,f) show J dxdy |E| 
and J dxdy |B| obtained from the cycle averaged fields E and B in the simulation box. They illustrate the magnitude 
of E and B along z in the bulk of the plasma. The plots show that close to the front surface and in the bulk of the 
plasma slab the magnetic field B is larger than the electric field E. At the rear of the slab E is larger than B. The 
electric field E at the rear of the slab points along the positive z-axis. At the front surface it points in the opposite 
direction. Furthermore, the relations E z 3> -Ex.y and B z <C B^ y are obtained in the simulation. Since the electric 
field is largest at the rear surface back surface acceleration of protons is most efficient in this simulation. 

Magnetic and current filaments are very small initially. There have been efforts to describe this early stage with 
the help of the Weibel theory |9pC[| . When two filaments approach each other the magnetic fields of the latter tend 
to cancel in the overlapping region leading to attractive forces between the filaments. This process is rapid in 3D and 
accompanied by a change of field topology. It cannot be recovered within the framework of the Weibel theory. 

Plots (a,b,c) of Fig. || show the cycle averaged current density filaments j z at different times for slices at z = 2.0 /xm. 
Plots (d,c,f) of the same figure show the total current and return current in z-direction obtained from I = J dxdy j z . 
The peak current density in the filaments shown in plots (a,b,c) is |j z | = 5.5 • I0 17 A/m 2 . The peak total current is 
about / = 7.0 • I0 5 A. Total current and return current are of equal magnitude while there are large discrepancies 
between both in individual filaments. The total current is largest close to the front surface of the slab. It declines 
in the bulk plasma. It has been split into its positive and negative components. To estimate the Alfven current 
Ia = 47reorn e c 3 /?7/e the quantities (3 = v/c and 7 = 1/ \J\ — (3 2 are required. The simulation yields (3 w 1.0 and 
7 w 2.0 for the central filament in plot (c). This gives I a ~ 34.0 kA. However, the total current in the central filament 
is about 200 kA which exceeds the Alfven current almost 6-fold. We note that the quantity (3 has been calculated 
from /3(x, t) = J d 3 pv z /(x, p, t)/c J d 3 pf(x, p, t) averaged over the cross-sectional area A of the central filament in 
plot (c). This yields (3 = J dxdy /3(x, £) / A. To make (3 as large as possible only the fastest electrons are included in 
the calculation with the requirement that the net current flowing in the filament is recovered from the latter. 

Figure || shows the formation of hollow, honeycombed density channels in the electron density. The same are 
obtained in the proton density (not shown in the present paper). Plots (a,b,c) of the figure show slices at z = 2.0 /xm 
for different times. Plot (d) gives a density slice at x = 2.0 fim. In plot (d) a central channel (blue color) surrounded 
by channels to the left and right of it is visible. These channels contain many hot electrons that exit the rear surface 
of the plasma slab. The energy density inside the channels is enhanced. The channels to the left and right of the 
central one coil around the latter. Hence, they are not fully visible in the plane at x = 2.0 /im. Comparison of Figs. 
|l] and Figs. || reveals that the hollow electron density channels are generated by the pressure of the magnetic field. 

In the context of current transport the singular current I s is of relevance. It is defined as the total current I 
minus the guiding center current I gc . The guiding center current comprises basically those electrons that gyrate 
around the magnetic field lines while the singular current contains mainly electrons that meander between field lines. 
The singular electrons are those that are capable of transporting large currents. The simulation yields approximate 
pressure balance between magnetic and thermal pressures on slowly varying time scales. Hence, we find 



diPij « (j x B),. , P iS = J d^ppiVjf . (1) 

The quantity Py in Eq. ([!]) is the pressure tensor which is calculated directly in the simulations. Rewriting Eq. (|l|) 
yields 



Equation (^J) is obtained from Eq. ([!]) assuming cylindrical geometry for a filament and rewriting Py with the help 
of the pressure Pj_ in radial direction normal to the z-axis. The magnetic field has been approximated by Bg and 
the current density by j z . Cylinder symmetry is approximately applicable to the situation found in the simulation. 
Hence, the calculation reveals the impact of geometry, pressure, and current density on the total current that can 
flow. Ambiguities inherent in the calculation of Ia as discussed above are avoided. The singular current I s is obtained 
with the help of Eq. (§) 

I s =I-I IJC = -^e c 2 ^-. (3) 

For a mono-energetic cold electron beam with electron density no and velocity v for which we have Pj_ = m e noV 2 
and j z = —eriQV the singular current I s in Eq. (||) yields the Alfven current Ia- These formulas have been obtained 
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in the context of Z-pinchcs Jll| . In hot laser plasma I s can become much larger than I a since hollow cylinder-like 
density channels form and a significant lateral pressure is present. As a consequence, larger total currents than the 
Alfvcn current arc possible in a single filament. From the simulation it is obtained that the singular current I s in the 
central filament of plot (b) in Fig. || is about 100 kA. Roughly this value is also obtained for / = J dxdy j z from the 
simulation by integrating over the cross sectional area of this filament. 

Figure [| shows the electron energy spectrum (a) , the integrated longitudinal and lateral charge and energy flows 
(b), and the electron number and energy densities in the plasma slab (c). The longitudinal flows have been split into 
positive and negative components. The return energy flow is negligible. For the lateral flows the integrated absolute 
values are taken. The fastest electrons obtained in the simulation have about 30 MeV while fractional absorption of 
the laser is about 35 %. The absorbed power of the laser is transported away by the energy flow q z shown in plot (b). 
The latter consists predominantly of fast electrons. The energy flow penetrates much deeper into the plasma than 
the charge flow but still decays rapidly. Plot (c) shows that the energy density is large only in the front layer of the 
plasma slab. Since the longitudinal and lateral charge flows are only sizeable where the magnetic field is large (see 
Fig. |l|) the implication is that the mutual interaction between current and magnetic field confines most of the total 
charge and energy flows to the front surface of the slab. Since time-averaged electric fields in the simulation are much 
smaller than time-averaged magnetic fields the transport regime is of magneto-hydrodynamic nature. 

In conclusion it has been shown that large charge flows generated by intense laser radiation decay into current 
filaments. The current filaments are isolated by magnetic filaments which reconnect until larger ones are obtained. 
This process proceeds until a dominating filament is obtained in the plasma. The magnetic pressure associated with 
the magnetic filaments leads to hollow density channels. Hollow density channels, however, can support total currents 
that exceed the Alfven current. This regime of transport is different from the one investigated by Bell et al. [ fl2| . 
Despite hollow channel transport severe charge and energy flow inhibition arc observed. We note that the Alfvcn 
current limit is difficult to apply in the context of a hot laser plasma since it has many ambiguities. We note that the 
plasma volume investigated is small due to numerical expense and that lateral boundaries are periodic. Hence, some 
care has to be taken when extrapolating the results reported here to FI relevant problems. 
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FIG. 1. Merging of magnetic filaments and cycle averaged E and B fields vs z in units Eq and Bo- Plots (a,b,c) show 
the planes z = 2.0 /im of the cycle averaged magnetic field B = (i? 2 + By + -B 2 ) ' 5 at different times. Plots (d,e,f) show 
E = J dxdy (El + E 2 + £ z 2 ) 5 (dashed) and B = J dxdy (B 2 + B 2 + B 2 ) 5 (solid) averaged over the lateral area and a laser 
cycle. The times are t = 71 fs (a,d), t = 125 fs (b,e), and t = 179 fs (c,f). The arrows in (a,b,c) indicate the direction of the 
cycle averaged magnetic field B. The parameters are Eo = 1.0 ■ 10 11 V/m and Bo = 8.76 • 10 2 Vs/m 2 . 



FIG. 2. Cycle averaged j z in units jo- Plots (a,b,c) show the planes z = 2.0 /im of the cycle averaged current density j z at 
different times. Plots (d,e,f) show I = Jdxdyj z split into flow (dashed-dotted) and return flow (solid) averaged over a laser 
cycle. The total current in the central filament (blue color) is approximately 200 kA. The times are t = 71 fs (a,d), t = 125 fs 
(b,e), and t = 179 fs (c,f). The arrows in (a,b,c) indicate the direction of the cycle averaged magnetic field B. The parameter 
is jo = 1.66 ■ 10 16 A/m 2 . 

FIG. 3. Hollow density channels. Cycle averaged electron density n e in units no. Plots (a,b,c) show the planes z — 2.0 /im 
at different times. Plot (d) shows the plane x = 2.0 /im. The times are t = 71 fs (a), t = 125 fs (b), and t = 179 fs (c,d). The 
parameter is no = 1.0 ■ 10 2S m -3 . The location of the hollow density channels coincides with magnetic filaments as is seen by 
comparing with plots (a,b,c) of Fig. |l| 

FIG. 4. Energy spectrum (a), total current and energy flows in units 7o and qo (b), total electron and electron en- 
ergy density in units no and eo (c). Plot (a) shows the electron energy spectrum obtained from all the particles in the 
plasma. The fastest electrons have 30 MeV. The thin solid lines of plot (b) show the integrated current and return current 
I = J dxdy j z . The thick solid line indicates I = J dxdy yj j| + jy ■ The thin dashed-dotted lines indicate the energy and 
return energy flows q z = J dxdy j d 3 pv z (cy m 2 c 2 + p 2 — mc 2 ) /. The thick dashed-dotted line shows the lateral energy flow 
q± — J dxdy J d 3 p y/v£~+v§ (cv/ m 2 c 2 + p 2 — mc 2 ) /. The solid line of plot (c) shows the integrated electron energy density 
e = j dxdy J d 3 p (c\J m 2 c 2 + p 2 — mc 2 ) /. The dashed-dotted line gives the electron number density n c — > J dxdyn c . Frac- 
tional absorption is about 35.0%. The time is t = 179 fs. The parameters are Io — 10 5 A, g = 10 12 W, eo = 10 4 J/m, and 
n = lO^nr 1 . 
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